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Abstract 

The impact of Hartree-Fock correlations on the nuclear momentum distribution is studied in a fully 
relativistic one boson exchange model. Hartree-Fock equations are exactly solved to first order in the 
coupling constants. The renormalization of the Dirac spinors in the medium is shown to affect the momentum 
distribution, as opposed to what happens in the non-relativistic case. The unitarity of the model is shown 
to be preserved by the present renormalization procedure. 
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It is well known that in a non-relativistic framework the momentum distribution of nuclear 
matter is not affected by the Hartree-Fock (HF) field. This arises because, due to general invariance 
principles pi, the non-relativistic self-energy cannot depend on spin in an infinite system: hence 
the single-nucleon wave functions are not modified and only the energy-momentum relation is 
affected by the medium. Of course correlations in the nuclear wave function beyond the mean 
field approximation are very important already at the non-relativistic level j|, ||, 0]. Due to such 
correlations, the momentum distribution is reduced for momenta below kp and the states above 
kp acquire small but finite occupation probabilities. 

Although the momentum distribution is not an observable, it is also true that over the years 
electron scattering reactions have frequently been expressed in terms of momentum densities. In 
recent work || we have evaluated the impact of mesonic correlations and meson-exchange currents 
(MEC) on the electroweak response functions within a fully relativistic, gauge invariant model. We 
have shown that the consistency of the theory necessarily implies the inclusion in the calculation of 
Hartree-Fock self-energy insertions. In order to deal properly with the divergencies associated with 
these diagrams, not only the energy but also the nucleon wave functions must be renormalized by 
the medium. As a consequence, in a relativistic HF framework, the momentum distribution is also 
modified for k < kp, since now the Dirac four-spinors describing the nucleons display new features 
and the self-energy becomes spin-dependent. The aim of this letter is to quantify this genuine 
relativistic effect in a one boson-exchange model for the NN interaction, while the corresponding 
observable consequences on the response functions were analyzed in depth in 

An unambiguous treatment of relativistic Hartree-Fock does not exist in the literature, since 
the presence of the Dirac sea requires (at least) the specification of a prescription to take it into 
account ||7j]. The approach we use is equivalent to that used in [[?], ||, where the nucleon proper 
self-energy is calculated in terms of positive-energy spinors only. This approximation is valid in the 
first iteration of a fully self-consistent calculation to which we confine ourselves in this work. This 
procedure was shown in @ to reproduce the non relativistic HF equations in the limit M — > oo, 
and it reduces the relativistic Hartree approximation to the Mean Field Theory. 

The proton momentum distribution of nuclear matter in the independent particle approximation 

is 

«(p) = E <»vMp)^ - k ) > (!) 

where kp is the Fermi momentum. Since we are focusing on symmetric nuclear matter the neutron 
and proton momentum distributions are equal. For a free relativistic Fermi gas, in momentum 
space the wave function describing a nucleon with momentum k and spin s is given by 



V>k>) = J v dv^ s {v)e'^ r = ^jL Ua (k,m) ^dre^-*** = ^ u s (p,m)5^ p (2) 

where V is the volume enclosing the system, u s (k, m) is the free Dirac spinor (fiu s = mu s ) and 
Efc = Vk 2 + m 2 is the free energy of the nucleon. We use the Bjorken and Drell ||] conventions for 
the spinor normalization, uu = 1. Therefore the wave function in coordinate space is normalized 
to one: f v drip^ s (r)ip kiS (r) = 1. 

The wave function of Eq. (Q), inserted in Eq. (|T|), yields the well-known result: 

n(p) = Vj2^-ul(P,m)u s ( P ,m)9(kp-p)5 pM = 2Ve(kp- P ) . (3) 
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In an interacting system, in the relativistic HF approximation, the above distribution is mod- 
ified, since the single-particle wave functions are renormalized by the interaction with the other 
nucleons in the medium. In this case the Dirac equation in the nuclear medium is given by 



[/>-m-£(P)]0 s (p) = O, 



(4) 



where <ft s (p) is the renormalized spinor and S(P) is the self-energy of a nucleon in nuclear matter. 
According to general symmetry properties £(P) can be written in the form [||, [l0| : 



£(P) = mA(P) + P(P) 7 oPo - C(P) 7 • P • 
Using the above decomposition the Dirac equation (Q) can be recast as 
[1 " C(P)] [ 7 o/o(P) - 7 • P - m(P)] MP) = , 

where the functions 

l-P(P) 



MP) 

m(P) 



1 - C(P) 
1 + A{P) 



Po 



m 



(5) 
(6) 

(7) 
(8) 



1 - C(P) 
have been introduced. 

Equation (jTJ) has the same structure as the free Dirac equation; hence for the positive-energy 
eigenvalue one has 

/ 2 (P)=p 2 + m 2 (P) , (9) 

which implicitly yields, using Eq. (|7j), the new dispersion relation for the renormalized energy 
Po = e (p) °f the nucleon in the nuclear medium: 

1 - C(P) 



Po 



p 2 + m 2 (P) 



1 - B{P) 

The corresponding positive-energy spinor reads (see Refs. []|, ||] for details) 



(10) 



4> s (p) = u s (p, m(p)) = JZ 2 (p 



E(p) +m(p) 
2m(p) 



1/2 



cr • p 



E(p) + m(p 



■Xs 



Z2(p)u s (p,m(p)) , 
(11) 

where the function m(p) of the three-momentum p is obtained from the Dirac mass in Eq. ^) by 
setting po = e(p): 

fn(p) = m(e(p),p) (12) 

and 

£(P) = /o(e(p), P) = \/p 2 + rh*(p) (13) 



represents the nucleon's Dirac energy. The field strength renormalization constant, \J Z 2 (p) , in 
Eq. (11) is obtained from the renormalized nucleon propagator || |ll]] and reads 



Z 2 (p) 



Res 



[1-C(P)][/ (P)-P(P)] po=e(p) 

dB m dA p 2 dC 
1 - B -p - m^-- h ^-t^— 

opo E op E opo 



(14) 



P0=e(p) 
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In the relativistic Hartree-Fock model the free spinors are used to compute the first approxi- 
mation to the self-energy. This is then inserted in the Dirac equation to get new spinors, and so 
on. This self-consistent procedure has to be dealt with numerically. 

The resulting momentum distribution is then obtained with the renormalized wave functions 



<MP) = / dr^ s (r)e-^ = 5,(p, m(p))«J k)P (15) 

Jv V e (p) 

and reads 



«(p) = E4, s (p)^ s (p)^(^-p) , (i6) 

where kp, m(p) and £?(p) are the nucleon's renormalized Fermi momentum, mass and energy, 
respectively. From Eqs. (flq. 15,11 ) the HF momentum distribution is then found to be 

n(p) = 2VZ 2 (p)9(k F -p) , (17) 

which clearly reproduces the free result in Eq. (||) for ^(p) = 1 and kp = kp. 

Note that the HF wave function with the spinor (ITl]) is not normalized to unity. Indeed 



/ *4 s (r)V k , s (r) = ^nt(k,m(p))n s (k,m(p)) = Z 2 (k) . (18) 
Jv ' E 

However the total number of nucleons must be conserved. This implies that the unitarity condition 

' /P f n(p) = 2V [ 7 ^-Z 2 (p)e(k F ~P) = Z (19) 



(2vr) 3 vtv J (2vr) 



must be fulfilled. Equation ( |l9| ) can be viewed as the procedure to fix the HF Fermi momentum 
hp, which can in principle be different from the free one. 

Here we consider the first-order correction to the momentum distribution arising from the HF 
series. We shall focus on mesonic correlations, induced by the exchange of n, p, u> and a, associated 
with the following interaction Lagrangian 12] 

C{x) = i>{x) { f^Vr • d^{x) + g p [ 7 M - ^^^] r ■ p^x) + g u -f^{x) + gM^)} 1>W . 

(20) 

Using this Lagrangian we compute the self-energy in OBE approximation. For each meson i = 
tt, p,u>,a, the corresponding self-energy functions Ai(P),Bi(P),Ci(P) are given in the Appendix. 
The total self-energy is obtained from Eq. (|) with A(P) = J2iM p ), B ( p ) = Ei-Bi(-P) and 
C(P) = J^jCj(P). While the pion and rho self-energies correspond to purely exchange (Fock) 
terms, the sigma and omega also have a direct (Hartree) contribution due to their isoscalar nature. 

The HF energy e(p), the solution of Eq. (|l0|), can be computed analytically to first order in the 
squared meson- nucleon coupling constant gf. For this purpose we note that the functions Aj(P), 
P>i{P) and Ci(P) are of order O(gf). Hence the following expansion of the Dirac mass in Eq. (B) 
holds: 

m(P) = m [1 + A(P) + C(P)] + 0{gf) . (21) 
Inserting this into Eq. ( ]lQj ) and expanding the righthand side to first order in gf , we get the equation 

po * E p + AE(P) , (22) 
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where 



AE{P) 



E r 



m 2 A{P) + E%B(P) - p 2 C(P) 



(23) 



Next we insert the value of po given by Eq. ( [22] ) inside the functions A(P), B(P), C(P) and expand 
them around the on-shell value po = E p , neglecting terms of second order in g 2 . We get 



A(P) ~ A(E p + AE, p) ~ A(Ep, p) = A (p) 



(24) 



and likewise B(P) ~ Bq(jp), C(P) ~ Co(p). Inserting these on-shell values into Eq. (|23| ) we finally 
obtain the HF energy to first order: 



E p + — m 2 A (p) + ^So(p) " P 2 C (p) = e(p) . 



(25) 



We proceed now by expanding as well the renormalized wave function, see Eqs. (15,11). For 

(26) 



this purpose we expand the Dirac mass in Eq. (21) around the on-shell energy 

m(p) ~m[l + A (p) + C (p)] 
and likewise the Dirac energy E(p) defined in Eq. ([I 

1 - B 



E(p) 



m 



l-C 



e(p)~Ep + — [A (p)+C (p)] 
E 



Moreover for the field-strength renormalization of Eq. (|l4|) we obtain 

Z 2 {p) ~l + a(p) 

with 

a(p) = B (p) + 



(27) 



(28) 



m 



2 dA 



dB p 2 dC 



Po=E p 



E p dp ' p dp Ep dp 
After some algebra the following first-order expression for the HF wave function is obtained 

A)(p) + C (p) £ P 70 - m 1 



(29) 



1 + m- 



E r 



2E r 



+ o«(P) 



u s (p,m)5^p . 



(30) 



The above expansion transparently displays the effect of the self-energy on the free wave func- 
tion. Indeed the second term in the square brackets of Eq. (30) corresponds to a negative-energy 
component with momentum p. Thus, within the OBE potential approach the renormalized HF 
spinors in the nuclear medium are characterized by two new elements with respect to the bare 
u s (p,m): the term (-Ep7o — m)u s (p, m), directly connected with the negative-energy components 
in the wave function, and the term a(p), arising from the field strength renormalization \J Z^pJ. 
However, the negative energy component does not contribute to the momentum distribution in first 
order, where one simply gets 



n(p) ~2V [l + a(p)]0(k F -p) 



(31) 
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The explicit expression for the first order expansion of the function a(p) = a(p) = J2i a i(p) 1S 
given in the Appendix. Note that the Hartree self-energy of the u and a does not contribute to a. 



The unitarity condition of Eq. (19) becomes 

2V J J^Z 2 (p)8(k F -p) = + ^ J k Jp 2 dpa{p) = Z , (32) 

which is certainly satisfied by kp = kp, because the function a exactly satisfies 

rkp 

dpp 2 a(p) = (33) 



o 

(see Appendix). We have numerically checked that there is no other value of kp for which the 
number of particles is Z. Therefore to first order in gf the Fermi momentum is not affected by the 
Hartree-Fock field. This means that the present calculation respects not only Lorentz covariance, 
but also unitarity. 

In Fig. 1 we plot the first order Dirac mass of Eq. (^) (top panel) and the first order momentum 
distribution of Eq. ([n]) (bottom panel) as functions of p/kp for kp=250 MeV/c. The separate 
contributions of the various mesons are displayed. In the present calculation we empirically account 
for the short-range physics through the meson- nucleon form factors -Fi(k) = A ' 2+k j which cut off 

the nucleon- nucleon interaction in a spatial region of size ~ 1/Aj . Actually, for sake of simplicity, 
we have approximated their effect by multiplying the self-energy associated with each meson by a 
constant factor (1 for the pion, 0.9 for the sigma, 0.5 for the omega and 0.4 for the rho): the form 
factors are indeed slowly varying functions of the meson momentum in the integration domain. 
The figure shows that the most sizable contribution to the Dirac mass arises from the cr-meson, 
which reduces the mass by about 30%, whereas the impact of the other mesons is at most 10% (in 
particular the pion induces a negligible increase of the mass) : the total effect in the present model 



is a reduction of the mass by about a factor 0.6, in accord with the findings of Refs. []7|, |10|, |l3j. It 
is also remarkable that the p-dependence of the Dirac mass is almost negligible. 

As far as the momentum distribution is concerned, it appears that the a, carrying an attractive 
interaction, induces a depletion of the baryonic density at low momenta and an enhacement of the 
latter in the vicinity of the Fermi surface, in contrast with the effect of the other mesons. It is 
interesting to note that the size of the Fock contribution in the momentum distribution decreases 
as the meson mass increases. This is in agreement with the fact that, at least for quasielastic 
inclusive electroweak responses modeled as we do here, the forces carried by the heavier mesons 
can be reasonably well approximated by four-fermion point interactions. In this case the HF 
approximation can be expressed as a linear combination of Hartree terms, which, as previously 
mentioned, do not affect the momentum distribution. Furthermore, and notably, the contributions 
arising from p, a and u cancel almost exactly. Thus the net effect of the full interaction coincides 
with the one obtained with the pion alone, and it amounts to an increase of the nucleon momentum 
density by about 1% for p ~ and to a decrease of it by almost the same amount for p ~ kp, 
in such a way that the number of nucleons is conserved, according to the unitarity condition of 
Eq. (JO]). It is worth noticing that the reduction of the momentum distribution near the Fermi 
surface due to relativistic HF correlations is of the same size as the one arising from short-range 



correlations of Jastrow type only, found in Refs. [14, 15 1 in the spectroscopic factors of quasi-hole 



valence states. Note, however, that this effect is very small compared with that expected from a 
more sophysticated non-relativistic modeling of short-range correlations 0], although one cannot 
make this statement with certainty in a relativistic context, since a relativistic version of Brueckner 
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Figure 1: The ratio m(p)/m of Eq. (|2q ) (upper panel) and the momentum distribution per unit 
volume of Eq. (lower panel) are plotted versus the nucleon momentum divided by the Fermi 
momentum (/cp=250 MeV/c). The results obtained by taking into account one single meson (pion: 
dashed; rho: double-dashed, omega: dotted; sigma: dot-dashed) are displayed together with the 

=770 MeV/c, m, 



total result (solid). The mesonic parameters are [24]: m n 
=550 MeV/c, ^/4vr 2 



=139.6 MeV/c, 



MeV/c, m a 
MeV/c, A p = 



13.6, ^/4vr 2 =0.84, a p =Q.l, ^/4vr 2 =20, ^/4tH=7.78, A w = 



=782 
1720 



1310 MeV/c, A w =1500 MeV/c, A a =2500 MeV/c . 



HF is even more challenging to carry out than relativistic HF and both constitute work for the 
future. 

In Fig. 2 the same observables displayed in Fig. 1, the contributions of all the mesons being 
included, are shown for three values of the Fermi momentum: kp = 200 (solid), 250 (dashed) and 
300 (dotted) MeV/c. It appears that the effect of the mesonic HF field on both the Dirac mass 
and the momentum distribution increases with the density and it is roughly proportional to kp. 
However the origin of this dependence is different in the two cases, due to the different role played 
by the various mesons. Indeed the Dirac mass essentially stems from the a and u, whereas the 
baryonic momentum density is significantly affected by all the four mesons. 

A simple analysis of the fcp-dependence of the OBE contributions to m(p) and n(jp) can be 
performed through an expansion in the small parameter rjp = fcp/m, whose typical value is ~ 1/4. 
Such an expansion has been successfully applied to the study of inclusive and exclusive electron 
scattering both for free and correlated nuclear systems ][l(], 17, 18, 19, 20]. Remarkably, the tjf 



expansion has been shown to be very useful for exploring the role of chiral pion dynamics in nuclear 



matter |2] 
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Figure 2: The ratio m(p)/m of Eq. (|2q ) (upper panel) and the momentum distribution per unit 
volume of Eq. (31) (lower panel) are plotted versus p/ttF for fcp=200 MeV/c (solid), kp=250 MeV/c 
(dashed) and k F =300 MeV/c (dotted). 



When performing this expansion one should pay attention to the fact that the pion is much 
lighter than the other mesons: this induces a different fcF-dependence for the pionic contributions, 
since kp/m^ cannot be treated as a small parameter. In fact it is easy to show that the heavy 
mesons' contributions to m(p) go as k% (the pion contribution is negligible). On the other hand in 
n the pionic effect grows as k F , while the heavy mesons contribute as k%. We recall that the a, u> 
and p almost cancel in the momentum distribution (see Fig. 1). 

It is of importance to notice that the physics of real nuclei roughly corresponds to the range 
200 < kp < 250 MeV/c: here our prediction for the Dirac mass, rh, is close to the one for the 
effective mass [||, ^3|. It is only for larger kp that the two quantities start to differ substantially. 

Finally, it is also interesting to note that the momentum distribution n(p) coincides with the free 
one, n(p) = 2, for a value of p/kF, which is independent of both the specific meson and the value of 
kp (see Figs. 1 and 2). This finding can again be interpreted on the basis of the above-mentioned 
expansion, which shows that the function a vanishes for p ~ ^/3/5&f. 

Before drawing our conclusions we would like to address the issue of the relevance of our 
findings on physical observables. In this regard we have shown in Refs. 0, |6| that the effect on 
the electromagnetic response functions including pionic correlations due to the modification of the 
momentum distribution is negligible (see in particluar Fig. 13 in §). One could not say this a 
priori and so an important conclusion of that work plus the deeper understanding presented in this 
letter is that at the level where relativity is dealt with in a consistent way such correlations appear 
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under typical circumstances to be perturbatively small. 

In summary, in this letter we have presented a relativistic analysis of the single-particle proper- 
ties of nuclear matter in HF approximation within a meson-exchange model. In particular, we have 
focused on the role played by the pion, rho, omega and sigma on the Dirac mass of the nucleon and 
on the momentum distribution. Whereas the momentum distribution is not affected by HF field in 
a non-relativistic framework, in the relativistic case it is slightly modified due to the renormaliza- 
tion of the spinors. In this work we have quantified this effect to first order in the coupling constant 
where the HF equations can be solved analytically. Using this solution we have demonstrated that 
the field strength renormalization function exactly satisfies unitarity at this order. 

Moreover, whereas for the Dirac mass, as is well known, the effect of HF mesonic correlations 
amounts to about 30-40% and mainly arises from the a and u mesons, we have shown that in 
the momentum distribution a cancellation among the heavier mesons occurs and the total result 
basically coincides with the pionic contribution, which amounts to a 1-3% effect, depending upon 
the density. 
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A Appendix 

The functions A, B, C of Eq. (|) can be expressed in terms of the integrals 

d 3 k 



I(P,m t ) 

L^(P,m t ) 

through the following relations 0: 

3tf 



(2tt> 
d 3 k 



-e(k F -k) 



(2vr) 3 
d 3 k 



6{k F - k) 
3 0(k F -k) 



1 


1 


2£ k (P 


- K) 2 - mf 


1 


K» 


2£ k (P 


- K) 2 - mf 


1 





MP) 
bap) 
c^p) 



P^(P,m n ) P 
I{P,rn^) - * 1 



2# k (P - K) 2 - mf ' 



2 - m 2 



2m 2 



M 

2 

M 

2 



I(P,m n ) 



I(P, m„ 



P^(P,m n ) , P 2 -m 2 L (P,m n ) 



+ 



m* 2m 2 po 

P^(P,m n ) P 2 -m 2 L 3 (P,m 7r ) 



+ 



2m 2 



for the pion, 



M p ) 



(2 + 3a p + 3oJ)/(P, m p ) - 3a p (l + a p ) m P ] 



(34) 
(35) 
(36) 

(37) 
(38) 
(39) 



J L is parallel to p since, choosing p along the z axis, the azimuthal integration in Eq. (B5h yields L\ — L2 — 0. 
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3a 2 

+ ^(P 2 -m 2 )I(P,m p ) 



B p {P) = 6g 2 ^3a p + 2a 2 p )I(P,m p ) 

2 



a 
m 



2 

P -P 

2 V 



2L»(P,m p 



L°»{P, 



Po 



1 + 3a + a 



«7 



L (P,m p 



C P (P) = 6g 2 J(3a p + 2a 2 )I(P,m p 



2 r > 



2L^(P,m p ) 



1 + 3 °P + °?-2^(^- m: 



L 3 (P,m p ) 
P 



for the rho, 





MP) 




B U (P) 




cup) 


and 






= -9l 


B«(P) 


2L 

= 9a 


CAP) 


= 9a 



kp L (P,m u 



3p -K 2 m 2 
2 L z (P,m UJ 
P 



Po 



r/n \ 1 /, _ 21 k F + E F \ 

IIP, m a ) H 5 — ^ fcir^ir — m In 

7r z m£. V m J 



Po 



P 

for the sigma. 

The corresponding expression for the functions «j (see Eq. (|29|)), with a = ^i=7r p< 



«i(p) 



U f kF ^k 2 E k -E, 



dk 



4tt 2 E p J E k7 2 (p,k)-4p 2 k 2 



fiip, k) 



where 



7i(p, fc) = (Ep - £ k ) 2 - p 2 - k 2 - m 2 = 2m 2 - m 2 - 2E p E k 
and the functions fi(p, k) are defined as 



Uip, k) 
f P (p, k) 



2 ^ ^ 2 m 2 p ,m 2 



2(1 + 6a p + 4a;) + at— £ + 4 



[2 + 6a p + a£-f 



p m 2 m 2 
2 m l\ 7 p (p> k) - 4p 2 k 2 Jp (p, k) + 2kp 



4kpm p 



IpijP, k) - 2kp 



UiPi k) 

fa(P, k) 



n , n m2 , lt{P,k) - 4p 2 k 2 ^(p, k) + 2kp 
1 + 2 — tt H — Ti m 



iri- 



2kpm 2 



1 _ 4^ + ^(P' fc ) ~ 4 P 2 k 2 ln >(P> *0 + 2k P 



mt, 



4kpm 2 



7 a (p, fc) - 2fcp 



The meson-nucleon form factors have been neglected for simplicity. Their impact on the results is 
discussed in the text. Using the above expressions, the unitarity condition of Eq. ( |33| ) follows since 
the functions fi(p,k) are all symmetrical under the exchange of p and k, hence 



k f 



dpp 2 ai(p) 



™%9i 
4tt 2 



dp 



kp p^k^ 
dk 



El, — Er 



E p E k Tf(p,fc)-4p2fc2 



fi(p,k)=0 



(55) 
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